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[t]  For  studies  of  heat  transfer  in  the  upper  ocean,  the  vertical  distribution  of  solar 
radiation  (ESR)  in  the  shortwave  domain  plays  an  important  role.  In  earlier  studies,  a  sum 
of  multiple  exponentials  was  used  to  describe  the  vertical  transmittance  of  ESr.  For  those 
exponential  terms,  an  attenuation  coefficient  for  each  term  is  assigned,  and  those 
attenuation  coefficients  are  assumed  to  be  vertically  constant.  Furthermore,  those 
attenuation  coefficients  are  empirically  modeled  as  functions  of  chlorophyll  concentration 
(Chi)  to  cope  with  varying  water  properties  of  the  oceans.  Since  attenuation  coefficients  of 
ocean  waters  are  generally  determined  by  components  more  than  Chi  alone,  we  developed 
a  generalized  model  bypassing  the  use  of  Chi  on  the  basis  of  extensive  numerical 
simulations.  In  this  new  model,  vertical  transmittance  of  ESR  is  separated  into  two 
exponential  terms,  with  one  for  the  contributions  of  wavelengths  below  700  nm  (the 
visible  domain,  £Vis)  and  one  for  wavelengths  above  700  nm  (the  infrared  domain,  Zsfr). 

An  attenuation  coefficient  is  assigned  for  each  of  the  two  exponential  terms.  Unlike  earlier 
models,  these  attenuation  coefficients  vaiy  with  depth  as  expected.  Furthermore,  the 
attenuation  coefficient  for  the  visible  domain  is  modeled  as  a  function  of  water’s 
absorption  and  backscattering  coefficients.  Since  absorption  and  backscattering 
coefficients  of  the  world  oceans  can  be  adequately  derived  from  observations  of  ocean 
color,  the  model  we  developed  provides  an  effective  and  adequate  means  to  describe  the 
three-dimensional  variation  of  both  £Vis  and  ESR  in  the  upper  layer  of  the  world  oceans. 

Citation:  Lee,  Z.,  K.  Du,  R.  Arnone,  S.  Liew,  and  B.  Penta  (2005),  Penetration  of  solar  radiation  in  the  Upper  ocean:  A  numerical 
model  for  oceanic  and  coastal  waters,  J.  Geophys.  Res.,  110,  C09019,  doi:10.1029/2004JC002780. 


1.  Introduction 

[2]  For  studies  of  heat  transfer  in  the  upper  ocean,  it  is 
important  to  reliably  provide  shortwave  solar  radiation  (Esr, 
see  Table  1  for  tenns  used  in  this  text)  in  the  upper  water 
column  [Zaneveld  et  al. ,  1981;  Lewis  et  al,  1990],  Knowl¬ 
edge  of  this  vertical  distribution  of  Esr  requires  information 
of  ESR  at  the  surface  and  how  ESR  attenuates  with  depth  (z). 
The  total  Es r  incident  on  the  water  surface  can  now  be 
estimated  with  models  or  measurements  made  from  space 
[Gregg  and  Carder,  1990;  Mueller  et  al.,  2004].  This  energy 
can  be  considered  as  two  separate  portions:  energy  with 
wavelengths  less  than  700  nm  (visible  domain,  A  vis),  and 
energy  with  wavelengths  greater  than  700  nm  (infrared 
domain,  EiR).  The  attenuation  coefficient  for  Em  is  almost 
completely  dominated  by  water  absorption;  thus  this  atten- 
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uation  can  be  viewed  as  invariable  to  changes  of  water 
constituents  [ Morel  and  Antoine,  1994].  Also,  because  of 
the  enormously  large  absorption  coefficients  of  water  mol¬ 
ecules,  this  portion  of  solar  radiation  is  quickly  lost  in  the 
upper  3  m  [ Morel  and  Antoine,  1994], 

[3]  £Vis  encompasses  the  wavelengths  shorter  than 
700  nm.  The  pioneer  study  of  Zaneveld  et  al.  [1981]  and 
subsequent  studies  [Chang  and  Dickey,  2004;  Kara  et  al., 
2005;  Kirk,  1988;  Lewis  et  al.,  1990;  McClain  et  al.,  1996; 
Morel  and  Antoine,  1994;  Ohlmann  et  al.,  2000;  Schneider 
and  Zhu,  1998;  Siegel  et  al.,  1995]  have  demonstrated  that 
the  penetration  of  Evis  in  the  upper  layer  of  the  ocean  plays 
an  important  role  in  heat  transfer  atid  sea  surface  tempera¬ 
ture.  Unlike  the  attenuation  coefficient  for  EiR,  however,  the 
attenuation  coefficient  for  Ryis  changes  horizontally  with 
constituents  in  the  water  [Jerlov,  1976;  Morel  and  Antoine, 
1994;  Zaneveld  et  al.,  1981],  Detailed  horizontal  informa¬ 
tion  regarding  the  attenuation  coefficient  of  £Vis  is  required 
to  determine  the  global  distribution  of  £vis(z)- 

[4]  In  the  past  two  decades,  multiple  exponential  terms 
[Paulson  and  Simpson,  1977;  Simpson  and  Dickey,  1981; 
Zaneveld  and  Spinrad,  1980]  were  generally  used  to 
describe  the  vertical  distribution  of  AVis  or  ESR  (more 
reviews  can  be  found  in  work  by  Morel  and  Antoine 
[1994]  and  Ohlmann  and  Siegel  [2000]).  For  these  expo¬ 
nential  terms,  an  attenuation  coefficient  (or  attenuation 
depth)  is  assigned  for  each  term.  These  attenuation  coef- 
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Table  I.  Terms  and  Definitions 


Definition 

Units 

«(X) 

absorption  coefficient  at  wavelength  X 

m~l 

*b(X) 

backscattering  coefficient  at  wavelength  X 

m'1 

Chi 

concentration  of  chlorophyll 

mg  m~3 

Ed\  z) 

spectral  downwelling  itradiance  at  depth  z 

W  m~2  nm"'1 

downwelling  solar  radiation  for  wavelengths  in  the  range  of 

700-2500  nm,  =  j  is^X,  z)d\ 

W  m-2 

^VTs(Z) 

downwelling  solar  radiation  for  wavelengths  in  the  range  of 

350-700  nm,  =  j  EJ\,  z)d\ 

W  m-2 

TsrU) 

J  350 

downwelling  solar  radiation  at  depth  2,  =  £vis(z)  +  E\\ d?) 

WnT2 

ATd(X) 

diffuse  attenuation  coefficient  for  £d(X,  z )  at  X 

m"1 

&ws(z) 

attenuation  coefficient  for  £Vis(z)  between  0  m  and  depth  2 

m~! 

kid?) 

attenuation  coefficient  for  EiR(z)  between  0  m  and  depth  z 

m'1 

7vis(z) 

vertical  transmittance  of  Evis(z),  =  £vis(r)/£vis(0  ) 

- 

r,  rU) 

vertical  transmittance  of  £ir(z),  =  £iR(z)/£JR(0“) 

- 

n z) 

vertical  transmittance  of  £sr(z),  =  Tsaffi/TsatO-) 

- 

Tv  is 

f?vis(0  )/£sr(0  ) 

- 

IOP 

inherent  optical  properties,  mainly  referring  to  the  absorption 
and  backscattering  coefficients  in  this  paper 

m~' 

Pd 

subsurface  irradiance  reflectance 

- 

Pf 

surface  Fresnel  reflectance 

- 

0,  is,,) 

solar  zenith  angle  above  (below)  the  surface 

_ deg_ 

ficients  are  kept  vertically  constant,  but  horizontally  vary 
with  Jerlov  [1976]  water  types.  Dividing  waters  into  differ¬ 
ent  types  is  coarse  anti  qualitative,  where  the  attenuation 
coefficient  associated  with  a  water  type  is  a  rough  estima¬ 
tion.  Furthermore,  variations  of  attenuation  coefficients 
within  a  water  type  are  ignored,  and  sharp  steps  between 
adjacent  water  types  result.  In  the  last  decade,  advancement 
in  remote  sensing  of  ocean  color  by  satellite  has  provided 
fine-scale  measurements  of  optical  and  biological  properties 
for  the  world  oceans.  Therefore  models  that  incorporate 
these  valuable  observations  better  fit  the  studies  of  solar 
heating  on  global  scales. 

[5]  At  present,  advanced  numerical  models  [ Gordon  et 
al.,  1975;  Liu  et  at.,  2002;  Mobley,  1995;  Morel  and  Gentili, 
1991]  have  been  developed  to  accurately  simulate  the 
subsurface  light  field  with  known  water  properties  and 
boundary  conditions.  These  models,  which  are  excellent 
for  studying  the  general  trend,  are  impractical  for  three- 
dimensional  global  models  because  of  their  computational 
inefficiency.  Simple  and  efficient  models,  but  without 
significant  loss  of  accuracy,  are  required. 

[6]  On  the  basis  of  the  strategy  of  case  1  and  case  2  water 
separation  [Morel  and  Prieur,  1977],  simple  and  explicit 
models  have  been  developed  to  incorporate  satellite-derived 
chlorophyll  concentrations  (Chi)  into  the  description  of  the 
horizontal  and  vertical  distributions  of  solar  radiation.  A 
model  for  the  spectral  diffuse  attenuation  coefficients 
(A/d(\))  has  been  specifically  developed  for  the  spectrally 
resolved  £Vis,  or  so-called  downwelling  spectral  irradiance 
(Ed(\)),  with  values  of  Kd(\)  described  as  empirical  func¬ 
tions  of  Chi  [Morel,  1988;  Morel  and  Maritorena,  2001], 
When  Chi  values  are  provided  via  satellite  observations 
of  ocean  color  [Carder  et  al.,  1999;  Gordon  and  Morel, 
1983]  (see  also  Ocean  color  algorithm  evaluation,  available 
at  http://oceancolor.gsfc.nasa.gov/REPROCESSING/ 
SeaWiFS/R3/),  Kd{\)  of  ocean  waters  can  be  evaluated  for 
desired  wavelengths.  When  this  information  and  surface 
Ed(\)  are  provided,  Ed(\,  z)  at  desired  depths  can  be  easily 
calculated  [Gordon  et  al,  1980],  as  can  Ey]S(z),  the  spectral 


integration  of  EJ\,  z)  (see  Table  1).  Note  that  for  the 
ultraviolet  to  visible  domain,  the  number  of  wavelengths 
can  be  more  than  100  for  fine  spectral  resolution. 

[7]  Although  this  wavelength-resolved  approach  might 
be  required  for  studies  of  phytoplankton  photosynthesis 
[Sathyendranath  et  al,  1989a],  it  is  found  “superfluous 
when  dealing  with  the  heating  rate”  [ Morel  and  Antoine, 
1994],  and  computationally  expensive  for  fine-scale  three- 
dimensional  modeling  of  the  global  oceans.  To  increase  the 
computation  efficiency,  without  significant  loss  of  accuracy, 
efforts  have  been  made  to  simplify  the  process  of  modeling 
the  vertical  transmittance  of  Ey |S.  For  this  purpose,  trans¬ 
mittance  of  A’vis  is  usually  simplified  to  a  few  exponential 
functions,  with  a  partition  factor  and  an  “effective  attenu¬ 
ation”  coefficient  assigned  for  each  term  [Morel  and 
Antoine,  1994;  Ohlmann  and  Siegel,  2000],  Furthermore, 
these  partition  factors  and  attenuation  coefficients  are 
expressed  as  functions  of  Chi,  and/or  cloudiness  ratios 
[Morel  and  Antoine,  1994;  Ohlmann  and  Siegel,  2000]. 
Therefore  when  Chi  values  and  other  auxiliary  information 
(such  as  Sun  angle)  are  available,  vertical  distribution  of 
E\ns  can  be  quickly  provided  and  efficiently  incorporated 
into  physical  models. 

[8]  However,  these  models  or  strategies  work  only  for 
case  1  waters,  where  all  optical  properties  are  determined  by 
Chi  alone  (with  solar  zenith  angle  explicitly  or  implicitly 
included)  [International  Ocean-Colour  Coordinating 
Group,  2000;  Morel,  1988].  Since  ocean  waters  have 
various  combinations  of  biological  and  chemical  constitu¬ 
ents  [Bricaud  et  al.,  1981;  Stramski  et  al.,  2001],  the 
globally  averaged  bio-optical  relationship  could  result  in 
big  uncertainties  for  a  specific  region  or  a  specific  time. 
Fundamentally  this  is  due,  first,  the  optical  properties  of 
ocean  waters  being  generally  determined  by  at  least  three 
independent  variables  [Sathyendranath  et  al.,  1989b]:  phy¬ 
toplankton,  colored  dissolved  organic  matter  (CDOM),  and 
suspended  sediments.  In  many  cases  knowledge  of  chloro¬ 
phyll  concentration  alone  is  inadequate  to  accurately  quan¬ 
tify  the  attenuation  coefficients,  an  optical  property.  Second, 
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even  for  chlorophyll  dominated  waters,  the  attenuation 
coefficient  can  differ  significantly  from  one  time  to  another 
or  one  place  to  another  because  of  the  “package  effect” 
[Bricaud  and  Morel,  1986;  Kirk,  1994],  These  natural 
variations  put  significant  limitations  on  the  attenuation 
coefficients  calculated  on  the  basis  of  Chi,  and  then  on 
the  resulting  vertical  distribution  of  Avis- 

[9]  To  avoid  such  limitations  associated  with  Chl-based 
models,  efforts  have  also  been  made  to  model  the  vertical 
transmittance  of  AVis  as  a  function  of  water’s  optical 
properties.  In  an  earlier  study  based  on  field  measurements, 
Zaneveld  et  al.  [1993]  modeled  the  attenuation  coefficient 
of  Av is  as  a  linear  function  of  the  diffuse  attenuation 
coefficient  at  490  nm.  Dividing  £Vis  of  the  upper  water 
column  into  three  consecutive  layers,  separate  empirical 
coefficients  were  developed  for  each  layer  [ Zaneveld  et  al., 
1993],  These  coefficients,  however,  remain  the  same  within 
each  layer,  which  thus  cannot  represent  the  steeper  than 
exponential  reduction  of  £Vis  in  the  subsurface  layer 
[Zaneveld  and  Spinrad,  1980].  Later,  Barnard  et  al. 
[1999],  on  the  basis  of  30  measurements  made  around 
California  waters,  modeled  the  attenuation  coefficient  of 
is  vis  as  a  polynomial  function  of  depth  and  absorption 
coefficient  at  490  nm.  This  empirical  model,  however,  did 
not  consider  the  effects  of  scattering  and  was  based  on  a 
relatively  narrow  range  of  water  properties.  Its  applicability 
to  broader  range  of  waters,  especially  where  scattering  by 
suspended  sediments  could  be  important  in  determining  the 
light  field,  is  unclear  [Barnard  et  al.,  1999]. 

[10]  Realizing  that  water’s  absorption  (a)  and  backscat- 
tering  coefficients  (bb)  can  be  routinely  measured  in  the 
field  and/or  well  retrieved  from  ocean  color  remote  sensing 
[Hoge  and  Lyon,  1996;  Lee  et  al.,  2002;  Loisel  et  al.,  2001; 
Lyon  et  al.,  2004],  it  would  be  very  useful  if  a  generalized 
model  for  the  attenuation  of  £Vis(z)  can  be  developed  that 
uses  values  of  a  and  6b  as  inputs.  To  reach  this  goal,  we 
carried  out  extensive  numerical  simulations  by  HydroLight 
[Mobley,  1995;  Mobley  et  al.,  2002],  and  developed  a 
simple  model  for  K v/s(z),  the  vertical  attenuation  coefficient 
of  Avis,  with  a,  bh,  and  z  as  variables.  Combined  with  the 
model  for  AiR(z),  a  model  with  two  exponentials  is  devel¬ 
oped  for  the  vertical  distribution  of  ASr(z).  With  this  model 
and  images  of  a  and  bb  provided  by  satellite  remote  sensing, 
three-dimensional  ASR(z)  in  the  upper  ocean  can  then  be 
efficiently  calculated  and  implemented  into  models  for  heat 
transfer  studies.  Separately,  Avis(z)  can  be  incorporated  into 
photosynthesis  models  [Behrenfeld  and  Falkowski,  1997]  to 
study  primary  productivity  in  the  oceans. 

2.  Vertical  Transmittance  of  E IR 

[11]  As  in  the  work  by  Morel  and  Antoine  [1994],  ASR(z) 
is  partitioned  into  two  portions:  £Vis(z)  (350-700  nm)  and 
AIR(z)  (700-2500  nm): 


not  fundamentally  different  from  that  of  Morel  and  Antoine 
[1994],  In  addition,  the  split  at  700  or  750  nm  makes  no 
difference  in  simulating  ASR(z),  as  long  as  it  is  by  the  same 
approach,  though  a  slightly  different  set  of  model  constants 
would  result. 

[12]  The  vertical  distribution  of  the  two  components  in 
equation  (1)  are  described  by  two  separate  exponential 
functions, 


EVIS{z)  =  Ey,s{  0-)e-k™Mz  (2a) 


Eir{z)=Eir{  0-)e~k^z.  (2b) 


Here  K  Vik.z)  and  K1R{z)  in  equation  (2)  are  not  the  values  at 
depth  z  (so-called  local  value),  but  a  spatial  average  between 
0  m  and  z.  Also,  K visiz)  and  KJR(z)  are  spectrally  averaged 
diffuse  attenuation  coefficients  that  are  weighted  by  the 
spectral  irradiance  in  the  visible  and  in  the  infrared, 
respectively.  Both  Kvls  and  K1R  vary  with  depth  and  Sun 
angle. 

[13]  Since  water’s  absorption  coefficients  are  extremely 
large  in  the  infrared  domain  [Morel  and  Antoine,  1994],  Km 
is  considered  independent  of  changes  in  water  constituents. 
Also,  as  most  of  the  solar  radiation  in  this  domain  is 
absorbed  within  the  layer  of  top  3  m,  Morel  and  Antoine 
[1994]  treated  KIR  independent  of  z.  Here  we  moved  one 
step  further  to  resolve  the  vertical  variation  of  Km.  For  this 
development,  we  first  calculated  K IR(z)  as  was  done  by 
Morel  and  Antoine  [1994], 


( 


p2500 

/  £,/(\,(TyX 

J  700 


A«(z)=-  In 

/  Ed(\,  0~)e-o,4x)*/cos(A) 

KJ  700 


(3) 


with  0W  the  subsurface  solar  zenith  angle.  aw(X)  is  the 
absorption  coefficient  of  water  in  the  range  of  700- 
2500  nm,  with  values  taken  from  Smith  and  Baker  [1981] 
and  Morel  and  Antoine  [1994].  Ad(X,  0_)  is  estimated  from 
MODTRAN  4  with  standard  marine  aerosols  (more 
common  for  the  this  study).  The  effect  of  aerosol  type  on 
the  calculated  Km  is  small  (~10%  for  z  =  0.5  m  and  ~7% 
for  z  =  2.0  m)  even  when  Ad(X,  0~)  estimated  using  strong 
continental  aerosols. 

[14]  From  equation  (3),  values  of  kIR(z)  in  the  layer  of 
0.01-3.0  m  (for  z  deeper  than  3  m,  Em  is  negligible) 
are  calculated,  with  Sun  angle  (0n)  varying  from  0-60° 
(equivalent  to  0W  from  0  to  40°)  from  zenith.  From  these 
Kir{z),  we  found  that  K IR(z)  can  be  described  as 


KIR{zA)  = 


0.560  + 


2.304 


(0.001  +zf 


(1  +0.0020,).  (4) 


Asr  (z)  —  Avis  (z)  +  Air  (z) .  ( 1 ) 

The  split  is  selected  at  700  nm  instead  of  750  nm  [Morel 
and  Antoine,  1994]  to  accommodate  outputs  and  models 
commonly  used  in  physical  and  biological  studies  [Gregg 
and  Carder,  1990;  Kirk,  1994;  Morel,  1978;  Zaneveld  etal, 
1993],  This  difference  is  simply  a  practical  choice,  and  is 


Figure  1  presents  examples  of  equation  (4)  modeled  KIR{z) 
versus  equation  (3)  calculated  KIR{z),  with  an  average  error 
of  ~3.0%.  Clearly,  for  a  change  of  2  orders  of  magnitude  of 
KIR(z)  in  the  surface  layer,  equation  (4)  represents  the 
vertical  variation  of  K1R(z)  very  well. 

[15]  Studies  [Pegau  et  al.,  1997;  Pegau  and  Zaneveld, 
1993]  have  found  that  values  of  aw(X)  in  the  infrared 
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inherent  optical  properties  (IOP)  [Preisendorfer,  1976],  A 
comprehensive  model  is  thus  needed  to  describe  the  change 
of  K vis  for  different  optical  properties  and  different  depths. 
To  develop  such  a  model,  concurrent  data  containing  wide 
range  of  KVI^z)  and  the  corresponding  optical  properties  are 
required.  In  an  approach  similar  to  Ohlmann  and  Siegel 
[2000],  we  used  HydroLight  to  generate  appropriate  data 
sets.  From  the  HydroLight  data  sets,  we  developed  a  model 
to  describe  K yjs  as  a  function  of  water’s  inherent  optical 
properties  and  depth.  Such  a  model,  combined  with  avail¬ 
able  optical  properties  from  satellite  measurements,  can  be 
applied  for  studies  of  heat  transfer  and  photosynthesis  in 
global  oceans. 

4.  Numerical  Simulations  to  Obtain  Kyis(z) 


Figure  1.  Example  of  modeled  K[R(z\  (equation  (4),  open 
diamonds)  compared  with  calculated  k m(z)  (equation  (3), 
solid  diamonds),  with  the  Sun  at  50°  from  zenith. 


domain  vary  with  temperature  and  salinity,  with  the 
temperature  effect  greater  than  the  salinity  effect.  This 
temperature  effect  is  found  generally  within  ~10%  of 
aw(X)  for  a  temperature  range  of  0°C  to  30°C,  and  the 
larger  effects  only  happened  at  a  few  wavelengths  (such  as 
750  nm)  [see  Pegau  et  a!.,  1997].  Since  Km(z)  represents 
the  spectrally  integrated  attenuation  in  the  wavelength  range 
of  700-2500  nm,  the  temperature  effect  at  a  few  wave¬ 
lengths  is  considered  insignificant  to  kiR{z).  Therefore 
values  presented  in  equation  (4)  should  be  applicable  to 
waters  with  temperature  and  salinity  widely  observed  in  the 
oceans. 

3.  Vertical  Transmittance  of  fvxs 

[i6]  For  the  visible  domain,  kns  varies  not  only  with 
depth  (the  vertical  direction),  but  also  with  water  constitu¬ 
ents  (the  horizontal  direction).  In  order  to  cope  with  this 
three-dimensional  change  of  kyjS,  Morel  and  Antoine  [1994] 
expressed  the  vertical  transmittance  of  EVjs  (represented  by 
Tvis)  as 


Tvis(z)  — 


Evisjz) 

Ens[Q~) 


2 


i=  1 


(5) 


Here  v,-  and  k,  are  parameters  for  7Vis  and  are  constants  for 
different  depths.  For  these  multiexponential  models,  the 
values  of  v,  and  kj  are  further  expressed  as  5-order 
polynomial  functions  of  log(Chl)  [Morel  and  Antoine, 
1994],  Therefore  given  Chi  and  z,  values  of  7ViS(z)  can  be 
quickly  evaluated. 

[17]  As  stated  earlier,  water’s  optical  properties  are  gen¬ 
erally  determined  by  constituents  not  limited  to  Chi,  but 
also  by  CDOM  and  suspended  sediments  [Sathyendranath 
et  al.,  1989b].  Furthermore,  Chi  values  derived  from  remote 
sensing  of  ocean  color  are  not  as  reliable  as  those  of 
absorption  and  backscattering  coefficients  [Lee  et  al., 
2002;  Lyon  et  al.,  2004],  To  get  more  accurate  vertical 
distribution  of  £yis>  an  efficient  and  reliable  approach  is  to 
estimate  KVIS  (an  optical  property)  directly  from  water’s 


[is]  As  in  earlier  studies  [Morel  and  Antoine,  1994; 
Ohlmann  and  Siegel,  2000],  k Vidk)  is  first  derived  from 
numerically  simulated  £vis(z)  that  corresponds  to  a  wide 
range  of  optical  properties.  As  described  by  Ohlmann  and 
Siegel  [2000],  simulation  of  Ed(z,  X)  is  carried  out  by 
HydroLight.  Unlike  the  simulations  of  Ohlmann  and  Siegel 
[2000]  where  water’s  lOPs  were  determined  by  Chi  only, 
IOPs  in  our  simulations  were  simulated  with  varying  Chi 
and  independently  varying  CDOM  and  suspended  sedi¬ 
ments,  as  described  by  Lee  et  al.  [2002]  and  International 
Ocean-Colour  Coordinating  Group/Ocean  Colour  Algorithm 
Group  (Model,  parameters,  and  approaches  used  to  generate 
wide  range  of  absorption  and  backscattering  spectra,  2003, 
available  at  http://www.ioccg.org/groups/OCAG_data. 
html).  Later,  K  vik2)  is  modeled  as  a  function  of  total 
absorption  and  backscattering  coefficients.  This  way,  the 
concentrations  of  Chi,  CDOM,  and  sediment  are  no  longer 
required,  as  long  as  accurate  total  absorption  and  backscat¬ 
tering  coefficients  are  provided.  In  natural  waters,  for  the 
same  optical  properties,  it  is  common  to  have  different 
concentrations  of  Chi,  CDOM,  and  sediments. 

[19]  For  numerical  simulations  by  HydroLight,  numerous 
descriptions  can  be  found  in  the  literature  [Berwald  et  al., 
1995;  Lee  et  al.,  1998;  Mobley  et  al.,  1993,  2002;  Morel  and 
Loisel,  1998].  The  following  summarizes  the  input  settings 
carried  out  in  this  study. 

[20]  1.  The  downwelling  irradiance  at  sea  surface  from 
the  Sun  and  sky  is  simulated  by  the  spectral  model  of  Gregg 
and  Carder  [1990].  A  wind  speed  of  5  m/s  is  assigned.  Sun 
angles  were  at  10°,  30°,  and  60°  from  zenith,  respectively, 
with  a  clear  sky. 

[21]  2.  To  determine  the  light  field  within  the  upper  water 
column,  HydroLight  requires  information  of  water’s  IOPs. 
For  the  attenuation  coefficient  of  downwelling  irradiance 
(the  focus  of  this  study),  the  important  IOPs  are  the  total 
absorption  (a)  and  backscattering  coefficients  ( bb ),  as  this 
attenuation  is  relatively  insensitive  to  the  scattering  phase 
function  of  particles  when  values  of  a  and  by,  are  fixed  [Lee 
et  al.,  2005a],  Chi  values  ranging  from  0.03  to  30.0  mg/m3 
(20  scales  with  a  roughly  50%  increasing  rate)  were  used  as 
one  free  parameter  in  order  to  provide  a  wide  range  of 
optical  properties.  Similar  to  earlier  numerical  studies 
[Gordon,  1989;  Kirk,  1991;  Morel  and  Loisel,  1998],  all 
IOPs  were  kept  vertically  constant  in  this  study. 

[22]  3.  Total  absorption  coefficients  were  expressed  as  a 
sum  of  the  contributions  of  pure  water,  phytoplankton 
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Figure  2.  Examples  of  Kyj^z)  for  different  water  proper¬ 
ties.  The  numbers  in  the  inset  are  values  of  a(490)  and 
/ib(490),  respectively.  Squares,  diamonds,  and  circles 
represent  Kyj^z)  from  HydroLight  simulations,  while 
dotted  lines  are  models  from  equation  (7)  (see  text  for 
details). 


pigments,  and  CDOM  (including  detritus).  Values  of  pure 
water  absorption  coefficient  were  taken  from  Pope  and  Fry 
[1997],  with  the  pigment  absorption  coefficient  determined 
by  Chi  values  [ Bricaud  et  al.,  1995].  Similar  to  Lee  et  al. 
[2002],  five  different  CDOM  absorption  coefficients  for 
each  Chi  value  were  randomly  assigned,  with  the  ratio  of 
CDOM  absorption  to  pigment  absorption  at  440  nm  in  a 
range  of  0.5  to  5.0  (resulted  in  a  range  of  0.0 17-3. 2  m_1  for 
a(440)). 

[23]  4.  Scattering  is  separated  into  molecular  and  nonmo- 
lecular  (collectively  called  particle)  components  [Gordon 
and  Morel,  1983;  Morel  and  Gentili,  1993;  Sathyendranath 
et  al.,  2001],  with  the  molecular  scattering  coefficient  taken 
from  Morel  [1974].  As  described  in  Lee  et  al.  [2002],  the 
particle  scattering  coefficient  at  550  nm  varied  with  Chi  in  a 
random  fashion.  The  phase  function  for  particle  scattering  is 
the  average  from  Petzold’s  measurements  [Mobley  et  al., 
2002].  With  these  different  combinations  of  Sun  angles  and 
water  properties,  a  total  of  300  (20  (Chi)  x  5  (CDOM  and 
particle)  x  3  (Sun  angles))  HydroLight  mns  were  carried 
out. 

[24]  5.  HydroLight  simulations  were  performed  for  \  in  the 
range  of  350-700  nm  with  a  10-nm  spectral  resolution.  For 
the  same  set  of  optical  properties  and  boundary  conditions, 
sensitivity  tests  with  HydroLight  run  at  a  5-nm  resolution 
showed  negligible  (<1%)  effects  on  7Vis(z)  (not  presented). 

[25]  6.  Five  geometric  depths  (excluding  0  m)  were 
selected  for  each  HydroLight  run,  with  depths  spread  within 
and  beyond  the  euphotic  zone  [Kirk,  1994].  Before  initiat¬ 
ing  HydroLight,  a  preliminary  estimation  was  made  to 
decide  the  deepest  depth  (z<j)  for  each  simulation.  This  zd 
is  corresponding  to  ~10“"J  of  Ed(\,  0_),  with  X  the 
wavelength  with  lowest  attenuation  coefficient  (the  deepest 
penetration  wavelength)  for  that  spectrum.  This  0  to  zd 
water  “slab”  is  then  sliced  into  five  consecutive  layers.  For 
the  five  layers  with  z,-  equal  to  the  position  of  the  bottom  of 
each  layer,  Eyls(zi)/Eyls(0~)  span  a  range  of  ~  0-80%.  In 


these  HydroLight  simulations,  corresponding  to  the  differ¬ 
ent  optical  properties,  the  first  (zi)  (and  last,  z5)  depths  of 
the  layers  ranged  from  1 .8  (147.4)  m  for  the  low  absorption 
case  to  0.2  (13.7)  m  for  the  high-absorption  case. 

[26]  7.  Bottom  reflectance  and  inelastic  scatterings  (such 
as  Raman  scattering)  are  excluded  as  they  are  beyond  the 
scope  of  this  study.  Inelastic  scattering  makes  negligible 
contributions  to  heat  transfer  and  photosynthesis  in  the 
water  column  [Morel  and  Gentili,  2004], 

[27]  Note  that  the  light  field  is  actually  determined  by  the 
absorption  and  scattering  properties  [Mobley  et  al.,  2002], 
the  bio-optical  models  and  the  Chi  values  used  in  the 
HydroLight  simulations  merely  provide  a  wide  range  of 
inputs  for  the  numerical  simulations  [Kirk,  1991], 

5.  Variation  of  K  yiS  and  Its  Modeling 

[28]  From  HydroLight  simulated  A’visfo),  K  yIS(z)  is  cal¬ 
culated  through 


Figure  2  presents  examples  of  K y/s(z)  for  a  few  cases.  Not 
surprisingly,  for  vertically  homogeneous  water,  K Vjs{z)  still 
changes  a  lot  (greater  than  a  factor  of  two)  in  the  subsurface 
layer.  Also,  K vis{z)  differs  significantly  for  varying  water 
properties.  All  these  variations  are  the  fundamental  reasons 
for  the  multiple  exponentials  used  to  cope  with  the  steep 
reduction  of  £Vis(z)  below  the  surface.  When  a  vertically 
constant  value  is  used  to  replace  ATWjS(z),  it  will  result  in 
overestimation  of  EVis  for  the  near-surface  layer  and 
underestimation  of  EyiS  for  deeper  depths.  To  accurately 
describe  the  vertical  distribution  of  £Vis  that  is  needed  for 
studies  of  heat  transfer  and  photosynthesis,  it  is  important  to 
know  and  model  the  vertical  variation  of  K  WlS(z). 

[29]  For  each  individual  K  yjs(z)  profile,  we  found  that  it 
could  be  well  reproduced  by  a  two-parameter  model, 

Km(IOP,  z)  =  K  ( IOP )  +  .  (7) 

Here  K\  and  K2  are  model  parameters  with  K\  for  the 
asymptotic  value  at  greater  depths  and  AX  more  important  to 
the  subsurface  K  VIS  value.  IOP  here  represents  combina¬ 
tions  of  different  inherent  optical  properties.  K{  and  K2  in 
equation  (7)  do  not  vary  with  depth  but  vary  with  IOP.  The 
dotted  lines  in  Figure  2  show  equation  (7)  modeled  Kyj^z) 
for  those  examples.  Figure  3  presents  the  result  of 
HydroLight  Kyjfe)  versus  equation  (7)-modeled  K yjs(z), 
for  HydroLight  simulations  with  the  Sun  at  30°  from  zenith. 
For  the  simulated  data  that  K Vis(z)  ranged  from  0.04- 
1.6  m"1,  the  average  error  between  FlydroLight  Kns(z)  and 
equation  (7)-modeled  KVJS{z)  is  2.2%  (maximum  error  is 
6.4%).  Such  results  clearly  demonstrate  that  equation  (7)  is 
adequate  to  describe  the  vertical  change  of  Kyj^z).  Thus  the 
vertical  profile  of  Tvls(z)  can  be  modeled  by  a  single 
exponential  function  with  two  parameters  (Kj  and  K%). 

[30]  For  this  single  exponential  term,  however,  K yjs(z) 
varies  with  both  z  and  IOP.  For  global  models  to  incorporate 
the  smooth  distributions  of  water  properties  obtained  from 
satellite  observation  of  water  color,  it  is  necessary  to  spell 
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Figure  3.  KVJS(z)  from  model  (equation  (7))  compared 
with  Kyjs(z)  from  HydroLight  (30°  solar  zenith  angle), 
indicating  that  K j//.dX)  can  be  well  described  by  equation  (7) 
with  two  parameters. 

out  how  Kx2  vary  with  water’s  optical  properties,  as  well  as 
with  solar  zenith  angle.  For  K  vis(z)  obtained  from  Hydro- 
Light  simulations  with  the  Sun  at  30°  from  zenith,  Figure  4 
shows  the  variations  of  Kx>2  versus  total  absorption  coef¬ 
ficients  at  440  nm  and  490  nm.  These  results  show  a 
positive  correlation  between  Kl2  and  a(440,  490),  respec¬ 
tively.  Since  Kx<2  also  varies  with  backscattering  coeffi¬ 
cients  (not  shown),  values  of  K]2  (and  then  K  ns(z))  cannot 
be  well  predicted  from  a (440,  490)  alone  (see  Figure  4). 
These  features,  however,  suggest  that  Kl2  can  be  modeled 
as  functions  of  both  total  absorption  and  backscattering 
coefficients. 

[31]  Since  Klt2  are  cumulative  effects  of  absorption  and 
scattering  in  the  visible  domain,  there  could  be  many  ways 
to  model  K I>2  with  different  combinations  of  inputs.  A 
simple  approach,  as  suggested  by  Zaneveld  et  al.  [1993] 
and  Barnard  et  al.  [1999],  is  to  use  absorption  and  back- 
scattering  coefficients  at  one  spectral  band,  where  data  from 
satellite  measurements  are  also  available,  to  model  K w». 
For  this  approach,  we  tried  to  model  ff,,2  using  the  values  at 


440  and  490  nm  as  inputs,  respectively,  and  found  using 
values  at  490  nm  provided  a  slightly  better  result.  Therefore 
empirical  models  for  Kx<2  are  obtained,  with 

{IOP)  =  Xo  +  X,  (o(490))°'5+  xA(490)  (8a) 


K2(IOP)  =  (o  +  <>(490)  +  »90).  (8b) 

Now,  the  values  of  xo,i,2  and  c0,i,2  are  independent  of  both 
depth  and  water  properties,  though  they  may  still  vary  with 
solar  zenith  angle. 

[32]  To  observe  the  variations  of  KX  2  for  different  solar 
zenith  angles,  values  of  KX  2  from  FlydroLight  simulations 
with  the  Sun  at  10°  and  60°  from  zenith  were  analyzed. 
These  results  suggested  that  Kx  only  varies  slightly  for 
different  solar  zenith  angles,  as  expected  from  the  asymp¬ 
totic  theory  [Zaneveld,  1989].  More  variations  with  Sun 
angles  were  found  for  AT2-the  parameter  important  for 
subsurface  solar  radiation.  On  the  basis  of  these  Kit2  values 
from  HydroLight  simulations,  equation  (8)  is  expanded  to 
include  Sun  angle  dependence, 

Kx(IOP,  0a)  =[xo  +  X,  (a(490))°’5+  xA(490)]  (1  +  oo  sin(0a)) 

(9a) 

K2(IOP,  0a)  =  [Co  +  >(490)  +  >,,(490)1(0,,  +  a2  cos(0o)). 

(9b) 

Note  that  0„  is  the  solar  zenith  angle  above  the  surface.  With 
this  model,  the  final  mathematic  function  for  7Vis(z)  became 

n7S{IOP,z,  0„)  =  e-kvB(toPfMz '  (10) 

By  fitting  7Vis  values  from  HydroLight  simulations  with  the 
combination  of  equations  (7),  (9),  and  (10),  values  of 
parameters  Xo,i,2,  ‘Xu, 2,  and  004,2  are  derived  (provided  in 
Table  2)  by  least  squares  curve  fitting,  as  done  in  earlier 
studies  to  derive  model  parameters  [ Morel  and  Antoine, 


Figure  4.  Values  of  (a)  Kx  and  (b)  K2,  which  are  model  parameters  for  kV!s(z),  versus  absorption 
coefficients  at  440  nm  and  490  nm. 
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Table  2.  Values  of  the  Model  Parameters  for  Kvid.z) 


Parameter 

Value 

Xo.  m“‘ 

-0.057 

Xi 

0.482 

X2 

4.221 

So,  m  1 

0.183 

^1 

0.702 

-2.567 

0(0 

0.090 

Otl 

1.465 

0(2 

-0.667 

1994;  Ohlmann  and  Siegel,  2000].  Figure  5  presents 
equation  (10)  modeled  Tvls(IOP,  z,  0„)  versus  Tyls(IOP,  z, 
0„)  determined  from  HydroLight  simulations.  For  those  Tvis 
values  (limiting  to  the  range  of  ~0.001  to  0.8),  bigger  errors 
happened  at  7Vis  <  0.003,  where  the  effects  of  EVis  on  heat 
transfer  and  photosynthesis  in  the  water  column  are  small. 
For  7Vis  >  0.003,  the  average  error  is  ~9%.  These  results 
indicate  that  the  simple,  optical  property-based,  model 
(equation  (10))  is  adequate  for  describing  the  vertical  profile 
of  EyjS(z)  for  different  waters. 

6.  Improvement  Over  Chi-Based  Models 

[33]  On  the  basis  of  equations  (1)  and  (2),  the  vertical 
transmittance  (=  E,SR(z)/£'SR(0~))  for  shortwave  solar  radi¬ 
ation  is 

T(IOP ,  z,  0O )  =  Fme-k™VOPfMz  +  (1  -  Fns)e-k'*{lM\  (11) 

with  Kyis  and  Km  provided  by  equations  (7)  and  (4), 
respectively.  Here  Fyis  is  the  ratio  of  £Vis(0_)/£sr(0~),  and 
its  value  is  around  0.424  [Mobley,  1994]  when  the  solar 
radiation  is  split  at  700  nm.  Equation  (11)  is  similar  to  the 
widely  used  model  of  Paulson  and  Simpson  [1977,  1981]. 
The  attenuation  coefficients  in  the  two  models,  however,  are 
fundamentally  different.  The  attenuation  coefficients  in 
equation  (11)  vary  with  depth  as  suggested  from  radiative 
transfer,  while  they  are  vertically  constant  in  work  by 
Paulson  and  Simpson  [1977]. 


Figure  5.  7VIS(z)  from  model  (equation  (10))  compared 
with  7ViS(z)  from  HydroLight  for  three  Sun  angles. 


[34]  In  order  to  test  the  performance  of  this  model  and 
earlier  models,  a  “standard”  T(z)  data  set  is  generated, 
which  is 

T(z)  =  FvisTvis{z)  +  (1  _  Fvis)Tir[z),  (12) 

with  7Vis(z)  created  by  HydroLight  simulations  and  T1R(z) 
on  the  basis  of  equation  (3).  The  HydroLight  simulations 
were  run  with  Chi  ranged  from  0.03-30  mg/m3,  but  with 
completely  different  absorption  and  scattering  coefficients 
as  compared  to  those  described  in  section  4.  The  depths  are 
fixed  at  1,  2,  5,  10,  20,  and  50  m,  that  is,  within  the 
commonly  observed  surface  mixed  layer.  Also,  the  Sun 
angle  is  set  at  45°  from  zenith.  Figure  6  shows  the  values  of 
Chi  and  a(490)  that  are  used  as  inputs  to  these  HydroLight 
simulations  and  to  the  T(z)  models.  The  dotted  line  in  Figure 
6  represents  case  1  a( 490)  for  those  Chi  values.  Clearly,  the 
data  set  follows  the  pattern  commonly  observed  in  the  field 
[Gordon  and  Morel,  1983;  Morel,  1988]. 

[35]  There  are  two  earlier  models  that  can  provide  T(z)  for 
given  Chi  and  z.  One  is  that  of  Morel  and  Antoine  [1994] 
(MA94  hereafter),  where  T(z)  is  described  as  a  sum  of  three 
exponentials;  another  is  that  of  Ohlmann  and  Siegel  [2000] 
(OS00  hereafter),  where  T(z)  is  evaluated  from  a  sum  of 
four  exponentials.  For  this  standard  data  set,  since  values  of 
Chi  (so  do  a(490)  and  />b(490))  are  known,  T(z)  can  be 
predicted  from  these  models.  These  T(z)  can  then  be 
compared  to  each  other,  and  to  the  standard  data  set. 

[36]  The  MA94  model  has  the  same  mathematical  for¬ 
mulation  as  equation  (12),  but  the  wavelength  range  is  from 
300  to  2750  nm  and  the  separation  of  visible  and  infrared  is 
at  750  nm.  Since  solar  irradiance  at  the  sea  surface  in  the 
range  of  300-350  nm  and  2500-2750  nm  are  only  about 
1%  of  the  total  irradiance,  we  applied  the  model  parameters 
of  Morel  and  Antoine  [1994]  to  the  standard  data  set  without 
any  change.  The  split  at  750  nm  makes  the  Fy iS  value 
different  from  the  Fyls  value  split  at  700  nm,  and  a  value  of 
0.56  [Morel  and  Antoine,  1994]  is  used  when  applied  to  the 
MA94  model. 


Figure  6.  Values  of  a(490)  and  Chi  used  in  section  6  for 
model  testing.  The  dotted  line  is  case  1  o(490)  for  those  Chi 
values,  derived  on  the  basis  of  the  Jfd(490)  model  of  Morel 
[1988],  and  u(490)  is  approximated  as  0.9  A"d(490). 
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Figure  7.  T(z)  from  models  compared  with  T(z)  from  the  standard  data  set,  where  Chi  varied  from 
0.03-30.0  mg/m3  and  z  varied  from  1  to  50  m:  (a)  MA94,  (b)  OS00,  (c)  the  model  of  this  study,  and 
(d)  comparison  of  T(z)  of  all  models  for  Chi  <  3.0  mg/m3.  Both  MA94  and  OS00  used  Chi  value  to 
predict  T(z),  while  model  of  this  study  used  a(490)  and  bh(490)  to  predict  T(z).  These  Chi,  <z(490),  and 
bb(490)  values  are  independent  of  those  discussed  in  section  4. 


[37]  The  OS00  model  (equation  (3)  of  Ohlmann  and 
Siegel  [2000])  is  developed  for  the  ratio  of  net  irradiance 
flux,  that  is, 


vtt.s  ESr{z)-Eu{z) 

T{z)  =  ~e^( W~ 


(13) 


where  EJz)  is  the  upwelling  irradiance  at  depth  z,  and  f  (z) 
is  modeled  as  a  sum  of  four  exponentials  similar  to 
equation  (5).  Compared  to  T(z)  defined  by  Morel  and 
Antoine  [1994]  and  in  this  study,  there  is 


n*)  = 


T{z) 

(1  —  Pf  ~  Pn)(l  —  Pd) 


(14) 


Here  p F  is  the  surface  Fresnel  reflectance  and  pD  is  the 
irradiance  reflectance  of  water.  For  a  Sun  angle  of  45°,  pF  is 
about  0.03  [Morel  and  Antoine,  1994].  The  p/>  value  is 
generally  less  than  ~0.01  [Morel  and  Antoine,  1994],  We 
also  noticed  that  the  model  for  f  (z)  [ Ohlmann  and  Siegel, 
2000]  was  developed  for  a  wavelength  range  of  250- 
2500  nm,  while  the  wavelength  range  in  this  study  is  350- 


2500  nm.  Since  solar  radiation  in  the  wavelength  range 
shorter  than  350  nm  is  less  than  4%  of  the  total  at  sea 
surface  [Mobley,  1994],  the  influence  from  this  wavelength 
mismatch  on  T(z)  is  very  small.  Therefore  T(z)  for  the  model 
of  Ohlmann  and  Siegel  [2000]  is  approximated  as  1 .057v(z) 
on  the  basis  of  equation  (14)  and  T{z )  is  calculated  as  in  the 
work  by  Ohlmann  and  Siegel  [2000]  for  given  Chi  and  z 
with  cloud-free  model  parameters. 

[38]  Figure  7  compares  model-determined  T(z)  values  to 
the  standard  T(z)  values  for  all  Chi,  with  Figure  7a  for  MA94, 
Figure  7b  for  OS00,  and  Figure  7c  for  the  lOP-based  model 
(equation  (11)).  Furthermore,  Figure  7d  compares  all  mod¬ 
eled  T(z)  values  for  Chi  <3.0  mg/m3,  which  is  usually  the 
upper  limit  for  oceanic  waters.  All  comparisons  and  analyses 
are  limited  to  T(z)  greater  than  0.1%  (a  total  number  of  474 
left  from  the  600  “standard”  values). 

[39]  Clearly,  for  all  water  properties  and  depths  studied, 
the  T(z)  values  from  the  two  Chl-based  models  can  differ 
significantly  from  the  standard  T(z )  values.  Also,  although 
both  MA94  and  OS00  used  the  same  Chi  values  to 
determine  J\z),  the  results  from  both  models  are  quite 
different.  The  R2  (square  of  correlation  coefficient)  values 
in  log  transferred  domain  are  0.79  (slope  of  1 . 1 8)  for  MA94 
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(Figure  7a),  and  0.72  (slope  of  1.12)  for  OSOO  (Figure  7b). 
On  average,  the  percentage  error  for  T(z)  in  log-log  scale 
are  ~60%  from  the  two  Chl-based  models,  with  MA94 
performing  slightly  better.  For  Chi  <  3.0  mg/m3  (Figure  7d), 
the  performances  of  the  two  models  are  not  significantly 
better.  From  these  Chl-based  models,  the  euphotic  zone 
depth  (where  7Vis(z)  =  1%  [Kirk,  1994])  can  be  quite 
different,  and  could  be  significantly  deeper  or  shallower 
than  the  actual  depth.  Note  that  the  approximations 
made  regarding  MA94  and  OSOO  will  only  cause  a  small 
percentage  difference  in  the  calculated  T(z)  compared  to 
their  original  models. 

[40]  The  T(z )  values  predicted  by  the  IOP-based  model 
(equation  (11))  matched  the  standard  data  sets  very  well  for 
all  waters  and  for  all  depths  (with  an  average  percentage 
error  of  ~6%  in  log-log  scale).  The  R2  value  is  0.99  (slope 
of  0.99)  for  data  in  Figure  7c,  also  in  log  transferred 
domain.  Compared  to  the  results  from  Chl-based  models, 
such  results  demonstrate  the  robustness  and  stability  of 
models  on  the  basis  of  inherent  optical  properties.  Conse¬ 
quently,  the  euphotic  zone  depth  can  be  estimated  more 
accurately.  Note  that  both  OSOO  and  the  IOP-based  model 
were  developed  with  data  from  PlydroLight  simulations, 
and  that  the  test  data  set  here  is  independent  of  the  data  used 
for  the  development  of  both  models. 

[41]  Apparently  the  results  presented  here  were  not  in 
favor  of  Chl-based  models.  It  could  be  argued  that  it  is  not  a 
fair  comparison  for  MA94  and  OSOO,  as  in  the  development 
of  those  models  Chi  values  were  the  only  variable  used  to 
determine  the  optical  properties  of  the  water.  The  test  data 
set  here,  however,  deliberately  disturbed  that  picture  (see 
Figure  6),  with  optical  properties  determined  by  Chi  and 
other  optically  active  constituents.  So,  a  mismatch  between 
T(z)  from  the  Chl-based  model  and  T(z)  of  the  test  data  set  is 
not  surprising.  This  mismatch  does  not  disprove  the  vali¬ 
dation  of  MA94  and  OSOO  for  water  environments  where 
optical  properties  quantitatively  follow  the  case  1  relation¬ 
ship.  More  realistically,  the  results  presented  here  should  be 
viewed  as  clear  evidence  that,  if  water’s  optical  properties 
do  not  follow  a  predefined  empirical  bio-optical  relationship 
(e.g.,  the  dotted  line  in  Figure  6),  the  predicted  vertical 
transmittance  of  £Vis  (and  £sr)  will  be  subject  to  big 
uncertainties.  In  addition,  since  optical  properties  of  natural 
waters  are  indeed  determined  by  more  than  the  value  of  Chi 
[Sathyendranath  et  al.,  1989b;  Stramski  et  al.,  2001],  this 
comparison  clearly  demonstrates  the  limitations  and  uncer¬ 
tainties  of  the  Chl-based  models,  as  suggested  by  Morel  and 
Antoine  [1994], 

[42]  For  the  same  wide  range  of  water  properties, 
however,  the  IOP-based  model  (equation  (11))  provided 
significantly  better  results  in  predicting  T(z)  for  the  upper 
ocean.  Since  remote  sensing  of  ocean  color  can  provide 
adequate  measurements  of  water’s  optical  properties  for 
global  oceans,  clearly  the  IOP-based  model  can  greatly 
improve  the  study  of  heat  transfer  and  photosynthesis  in 
the  upper  layer  of  the  global  oceans.  However,  it  is 
necessary  to  keep  in  mind  that  so  far  this  newly  developed 
model  and  the  models  of  Morel  and  Antoine  [1994]  and 
Ohlmann  and  Siegel  [2000]  are  only  tested  with  numerically 
simulated  data.  For  validation  and  wide  applications,  it  is 
essential  to  test  the  models  with  high-quality  data  obtained 
from  field  measurements. 


[43]  The  above  compared  T(z)  values  are  for  a  few 
discrete  depths.  To  reveal  the  vertical  details  of  T(z)  in  the 
subsurface  layer,  Figure  8  presents  examples  of  T(z)  of  top 
10  m  for  four  commonly  observed  Chi  values  in  the  field. 
Figure  8  also  shows  T(z)  from  the  standard  data  set.  Clearly, 
since  Morel  and  Antoine  [1994]  did  not  try  to  resolve  the 
details  of  T(z)  within  the  top  1  m,  larger  errors  are  found  for 
near-surface  T(z)  in  the  MA94  model,  as  pointed  out  by 
Morel  and  Antoine  [1994].  The  models  of  Ohlmann  and 
Siegel  [2000]  and  this  study  tried  to  resolve  the  details  of 
near-surface  T(z).  In  addition,  since  the  attenuation  in  the 
infrared  domain  is  dominated  by  that  of  pure  water,  the 
results  between  T(z)  from  OSOO  and  T(z)  from  this  study 
show  good  agreement  for  the  near-surface  layer.  Note  that 
for  z  <  1  m,  changes  of  water  constituents  have  only  minor 
effects  on  T{z). 

7.  Discussion 

[44]  The  T(z)  values  from  the  IOP-based  model  used  only 
absorption  and  backscattering  coefficients  at  490  nm;  opti¬ 
cal  properties  at  other  wavelengths  were  ignored.  The  robust 
results  from  the  IOP-based  model  (see  Figures  7  and  8) 
suggest  that  knowing  the  absorption  and  backscattering 
coefficient  at  490  nm  is  adequate  for  estimation  of  the 
vertical  transmittance  of  £Vis  and  Esn,  as  implied  in  earlier 
studies  [Barnard  et  al.,  1999;  Zaneveld  et  al,  1993].  This  is 
somewhat  surprising,  because  for  studies  such  as  the  esti¬ 
mation  of  primary  production  detailed  spectral  information 
is  important  [Sathyendranath  et  al.,  1989a].  The  reason  for 
this  relaxed  requirement  in  spectral  information  is  that  the 
focus  of  this  study  is  the  spectrally  integrated  solar  radiation 
that  is  available  for  water  heating.  There  could  be  quite 
different  absorption  coefficients  at  350  or  410  nm  for  the 
same  absorption  coefficient  at  490  nm  [Lee  et  al.,  2005b], 
and  these  differences  will  result  in  different  irradiance  for 
those  wavelengths  at  a  specific  depth.  These  irradiances, 
however,  are  only  a  small  portion  of  the  total  solar  radiation, 
and  the  decrease  of  irradiance  at  one  wavelength  may  be 
compensated  by  an  increase  at  another  wavelength  at  that 
depth,  so  the  sum  of  the  irradiance  of  all  wavelengths 
(which  gives  £Vis  and  £sr)  may  not  vary  much.  Such  a 
scenario  is  not  true  for  depths  below  the  euphotic  zone 
where  irradiance  is  eventually  concentrated  to  a  narrow 
band;  but  the  contribution  of  solar  radiation  at  such  a 
spectral  band  and  at  such  depths  is  small  compared  to  the 
absorbed  solar  energy  in  the  upper  water  column. 

[45]  The  results  presented  here  show  significant  improve¬ 
ments  in  modeling  of  the  vertical  distribution  of  £Vis  and 
ESR  in  the  upper  ocean  by  utilization  of  IOPs.  Such  an 
approach  is  adequate  for  satellite  remote  sensing  of  ocean 
color  and  in  environments  where  water’s  optical  properties 
are  well  measured.  However,  if  Chi  value  is  the  only 
available  parameter  for  a  water  environment,  since  there 
are  wide  uncertainties  in  estimating  its  IOPs,  we  expect  no 
significant  difference  in  estimated  T(z)  from  either  IOP- 
based  model  or  MA94  or  OSOO.  It  is  recognized  that  the 
accuracy  of  IOPs  (a(490),  for  instance)  derived  from  ocean 
color  to  the  best  can  be  within  ~10%  [Lee  et  al.,  2002]. 
This  error  will  be  propagated^  to  the  calculated  Kns{z)  and 
then  Tns(z).  However,  since  K yisiz)  is  a  spectrally  averaged 
value  with  a  big  portion  determined  by  the  large  absorption 
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Figure  8.  Examples  of  jT(z)  for  the  top  10  m,  with  T(z)  calculated  by  three  models  and  from  the  standard 
data  set.  Values  of  Chl,  a(490),  and  Ab(490)  are  provided  in  the  insets. 


coefficients  of  aw(\)  in  the  red  wavelengths,  only  a  fraction 
of  errors  in  IOPs  will  be  transferred  to  Kv/S(z). 

[46]  It  is  necessary  to  point  out  that  the  model  developed 
here  has  primarily  demonstrated  the  strategy,  concept,  and 
advantages  of  modeling  the  vertical  transmittance  of  £v is 
and  ESr  using  inherent  optical  properties.  The  work  here 
presented  a  model  with  the  major  variables  such  as  IOPs 
and  Sun  angle  incorporated,  and  the  model  is  for  clear-sky 
conditions  with  marine  aerosols  and  a  well-mixed  surface 
layer  [ Mueller  and  Lange,  1989].  When  the  sky  is  cloudy, 
reduction  of  the  total  downwelling  irradiance  and  the  ratio 
of  the  visible  portion  (Fv jS)  need  to  be  adjusted  accordingly 
[Behrenfeld  and  Falkowski,  1997;  Ohlmann  and  Siegel, 
2000].  However,  a  recent  study  [Bartlett  et  al,  1998] 
pointed  out  that  clouds  have  only  a  small  spectral  effect 
on  surface  £Vis>  which  suggests  that  model  parameters 
developed  in  this  study  could  be  applicable  to  cloudy 
conditions,  though  a  slight  adjustment  for  9a  (or  subsurface 
solar  zenith  angle  as  in  work  by  Sathyendranath  and  Platt 
[1988])  is  necessary.  Similarly,  when  the  subsurface  light 
field  is  more  diffuse  as  a  result  of  capillary  waves  [Zaneveld 
et  al.,  2001],  an  effective  0a  might  also  be  required. 
However,  the  effects  of  slightly  incorrect  on  the  attenu¬ 
ation  of  Ayis  are  much  smaller  than  that  of  incorrect 
absorption  and  backscattering  coefficients. 

[47]  One  should  take  caution  when  applying  the  newly 
developed  model  to  depths  that  are  below  the  well-mixed 


surface  layer,  especially  if  the  optica]  properties  at  those 
depths  are  significantly  different  from  those  in  the  mixed 
layer.  For  such  conditions,  the  evaluation  of  downwelling 
irradiance  at  those  depths  will  be  changed,  and  the  proper¬ 
ties  retrieved  from  observed  ocean  color  might  also  be 
affected  [ Gordon  and  Clark,  1980;  Sathyendranath  and 
Platt,  1989;  Zaneveld  et  al.,  1998].  A  critical  hurdle  in 
remote  sensing  is  the  lack  of  instantaneous  knowledge  of 
the  vertical  distribution  of  the  optical  properties  when  a 
satellite  sensor  makes  measurements.  Historical  data  may 
provide  a  general  guidance  about  the  vertical  distributions, 
but  those  measurements  are  insufficient  and  the  vertical 
distributions  vary  spatially  and  temporally.  A  method  to 
adequately  provide  the  details  of  the  vertical  distribution  of 
the  optical  or  biogeochemical  properties  of  a  water  body  is 
needed.  Combining  measurements  from  satellite  sensors 
with  outputs  from  dynamic  oceanography  models  [Chai  et 
al.,  2002;  McGillicuddy  et  al.,  1995]  may  provide  useful 
alternatives  to  overcome  the  limitations  of  significantly 
stratified  upper  water  column. 

8.  Conclusions 

[48]  In  this  study,  an  innovative  model  is  developed  for 
describing  the  vertical  transmittance  of  solar  radiation  in  the 
upper  layer  of  the  oceans.  Similar  to  Morel  and  Antoine 
[1994],  our  study  explicitly  partitioned  the  solar  radiation 
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into  two  portions:  one  for  the  visible  domain  (EVisj  350- 
700  nm),  and  the  other  for  the  infrared  domain  (Air,  700- 
2500  nm).  For  the  two  portions,  two  exponential  functions 
are  used  separately  to  describe  their  vertical  transmittance. 
For  each  exponential  term,  an  attenuation  coefficient  is 
assigned.  For  Em,  its  attenuation  coefficient  ( KIR{z ))  is 
modeled  as  a  function  of  z  and  solar  zenith  angle.  For  Avis, 
its  attenuation  coefficient  (K  y/s(z))  is  modeled  as  a  function 
of  depth,  solar  zenith  angle,  and  water’s  optical  properties, 
after  extensive  numerical  simulations  using  HydroLight. 

[49]  Compared  with  earlier  Chl-based  models,  the  model 
developed  here  has  fewer  empirical  parameters  and  provides 
robust  and  stable  results  in  describing  the  vertical  distribu¬ 
tion  of  the  visible  portion  or  the  full  spectrum  of  solar 
radiation  in  the  subsurface  layer.  Since  the  model  requires 
«(490)  and  hb(490)  as  inputs,  accurate  retrieval  of  these  two 
optical  properties  from  satellite  observations  of  ocean  color 
is  required.  When  these  properties  are  available  for  world 
oceans,  the  model  can  be  used  to  provide  efficient  evalua¬ 
tion  of  AVis(z)  and  ASR(z)  in  the  well  mixed  upper  water 
column,  which  in  turn  can  be  incorporated  into  physical 
models  to  study  heat  transfer  [e.g.,  Kara  et  al,  2005],  Also, 
after  consideration  of  the  average  cosine  [ Sathyendranath 
and  Platt,  1988],  AVIS(z)  can  be  easily  converted  to  photo¬ 
synthetic  available  radiation  (PAR)  and  incorporated  into 
models  for  estimating  primary  production  of  the  water 
column.  In  addition,  the  K ns(z)  model  provides  an  easy 
and  reliable  tool  to  predict  the  light  level  at  desired  depths 
from  satellite  IOPs,  needed  to  plan  the  C14  incubation  for  in 
situ  measurements  of  primary  production  [Barnard  et  al., 
1999]. 
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